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Abstract 

Classical,  one-dimensional  theory  of  hydrodynamic  penetration  is  used  as  the  basis  of  establishing  simplified 
analytical  relationships  describing  energy,  momentum,  and  power  deposition  during  hypervelocity  impact  events. 
A  concise  overview  of  the  1-D  model  is  given  followed  by  a  select  grouping  of  terms  into  relationships  that  offer 
first-order  criteria  for  making  engineering  design  considerations  on  relevant  applications  and  assist  in  the  analysis 
of  experimental  observations.  Momentum,  energy,  and  power  deposition  are  found  to  be  proportional  to  second, 
third  and  fourth  power  exponents,  respectively.  These  analytical  terms  are  presented  for  constant  velocity 
gradient,  i.e.  fixed  length,  rods  as  well  as  linear  velocity  gradient  rods,  such  as  shaped  charge  jets. 

The  role  of  penetrator-to-target  density  ratio  is  then  examined  in  terms  of  the  backflow,  or  reverse  flow  of  1- 
D  penetration.  Again,  the  non-dimensional  ratio  of  penetrator-to-target  mass  density  is  used  to  compare  the 
relative  velocity  of  material  flow  during  penetration.  The  relationship  highlights  the  role  of  penetrator  materials 
for  achieving  desired  effects  in  these  hypervelocity,  terminal  ballistics  events. 

Albeit  the  relationships  are  derived  on  the  assumptions  for  hydrodynamic  processes,  their  generality  of  form 
and  ease  of  implementation  make  them  a  useful  first-order  description  for  engineering  insight  and  application  over 
a  broad  range  of  velocities. 

Keywords'.  Hypervelocity  Penetration,  Energy  Deposition  of  Jets,  Reverse  Flow,  Shaped  Charge  Warhead 


1.  INTRODUCTION 

The  Bernoulli  equation  has  long  been  used  to  produce  simplified,  but  insightful  analytical 
relationships  for  such  useful  terms  as  the  penetration  depth,  penetrator  erosion  rate,  and  interface 
velocity  of  the  penetrator-target  for  hypervelocity  and  shaped  charge  jet  impacts  [1,2].  Further  efforts 
[3,4]  expanded  on  the  fundamental  equations  to  establish  relationships  of  energy  and  momentum  flux 
and  deposition  rates  for  constant  velocity  and  linear  velocity  gradient  rods.  The  work  of  Foster  and 
Wilson  serves  as  the  foundation  of  the  derivations  within  this  paper.  A  brief  review  of  the  basic 
equations  is  necessary  for  completeness  and  establishing  variations  of  the  forms  for  further  study  and 
discussion. 
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1.1  The  Case  of  a  Constant  Velocity  Rod 


Derivations  are  first  made  for  a  constant  velocity  rod.  The  Bernoulli  equation,  with  its  assumption 
of  incompressible  and  inviscid  conditions,  is  applied  along  the  centerline  of  a  hypervelocity  projectile 
impacting  a  target,  Figure  1. 
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V  =  penetrator  velocity 
u  =  interface  velocity 


Fig.  1.  Penetration  terminology. 


The  steady-state  assumption  is  most  easily  understood  in  the  reference  frame  being  fixed  to  the 
penetrator/target  interface,  giving  the  streamline  energy  equality  and  solution,  Equation  (1). 


\pP{V-u)  =\pTu 2  having  interface  velocity 


V 


u  = 


1  +  y 


(1) 


With  term  y  defined  for  the  square  root  of  the  ratio  of  target  density  to  penetrator  density.  The  rate  of 
erosion  of  the  rod  increment,  dl/dt,  and  the  rate  of  penetration  (  advancement  of  the  interface),  dp/dt,  are, 


dl 

dt 


V-u 


dl  _  Vy 
dt  1  +  y 


(2) 


Using  the  differential  chain  rule,  the  increment  in  penetration  per  unit  rod  length  consumed  is  also 
found  in  terms  of  y,  knowing  that  the  interface  velocity  is  also  u  =  dp/dt , 


dp  dt 
dt  dl 


f  V  } 

"l  +  / 

U  +  r) 

l  Vy  ) 

dp  1 
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The  momentum,  p  ,  and  energy,  E,  terms  are  defined  for  the  constant  velocity  rod  having  a  fixed 
cross-sectional  area,  A. 


p  =  ppAlV  and  E  =  j  ppAlV° 


(4) 


Again,  the  chain  rule  is  used  to  find  the  momentum  and  energy  deposited  per  unit  penetration  path 
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length,  as  well  as  the  power  density  of  the  impacting  rod. 
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1.2  The  Case  of  a  Linear  Velocity  Rod 

The  derivation  of  penetration  relationships  for  the  linear  velocity  gradient  penetrator,  described  in 
Figure  2,  is  now  given.  This  domain  is  suitable  for  idealized  penetration  descriptions  of  shaped  charge 
jets  solved  in  its  virtual  coordinate  frame  [5]. 

The  velocity  of  the  penetration  interface  at  any  time  after  impact  is,  u  =  dx/dt ,  with  the  position  of 
the  penetrator  being  generally  described  by  the  ever-changing,  eroding  tip  velocity  in  the  target  as, 

x  =  Vt  (6) 


A  non-dimensional  position  response  is  obtained  by  integrating  the  penetration  interface  relation 
with  substitutions  of  Eq.  (6)  and  recalling  that  u  =  dp /dt  =  V /(l  +  y)  to  get, 


Xr, 


\{0J 


1 +r 


(V) 


Through  judicious  manipulation  of  Eq.  (6)  and  substitution  of  Eq.  (7),  the  velocity  is  obtained  in 
terms  of  initial,  known  values  and  virtual  origin  coordinates, 


y  _  f_  _  X  X0  h 

t  XQ  t  Q  t 


V  =  v„ 


f  \-r 

x 


(9) 


Vxo  J 


Elere,  the  area  is  no  longer  a  constant,  but  is  a  time-dependent  parameter  changing  with  length  as 
governed  by  the  velocity  gradient  between  the  leading  tip  and  the  tail  mass  points. 


A(t)  = 


^0 

pP(K,P  -vtail)t 


(10) 
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[penetration  path  coordinates] 
Fig.  2.  Linear  velocity  gradient  rod  penetration. 


Time  can  be  eliminated  from  the  area  relationship  through  the  initial  conditions, 
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Giving  area  as  a  function  of  spatial  dependency,  known  quantities,  and  virtual  coordinates, 

A(x)  = 


™o  VtiP 


f  y-U+r) 
X 


Pn(Kip  ~VtJX0 


Vxo  J 


(12) 


Using  the  relationships  derived  for  the  constant  velocity  case,  Eq.  (5),  and  substituting  for  area,  Eq. 
(12)  and  velocity,  Eq.  (9),  the  momentum  and  energy  relationships  deposited  per  unit  penetration 
distance  for  the  case  of  a  linear  velocity  gradient  penetrator  are  found  as, 
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And  the  power  input  function  is  then, 
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dE  _  m0V£ 
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'  x  'l 

-(i+4r) 

dt  W,ip  -Vtau)x0 

ll  +  r) 

VX0  ) 

(15) 


2.  DISCUSSION  OF  THE  1-D  RELATIONSHIPS 

The  momentum  and  energy  relationships  for  both  the  constant  and  linear  gradient  velocity  rods  are 
similar  enough  in  form  to  discuss  interchangeably. 

2.1  Constant  Velocity  Rod 

The  simpler  case  of  the  constant  velocity  rod  is  discussed  first.  Even  though  the  power  density  is 
mainly  dependent  (third-power)  on  velocity,  there  is  a  non-linear  contribution  from  penetrator  material 
density.  The  plot  of  Figure  3  extracts  the  velocity-dependence  and  examines  just  the  ppf/{  1+y)  term 
versus  material  density  for  rods  penetrating  a  target  of  density  f>r  =  7850  kg/m3,  e.g.  steel.  The  benefit 
of  this  input  energy  must  be  considered  along  with  transmitted  shock  strength  and  efficiency  as  coupled 
by  the  Hugoniot  match  for  the  penetrator/target  combination.  There  are  target  failure  mechanisms, 
albeit  not  covered  in  such  a  simplistic  1-D  theory,  which  might  be  activated  using  a  rapid  deposition  of 
energy.  Mechanisms  of  brittle  fracture,  pore  collapse,  spallation,  phase  change,  and/or  crush-up  are  but 
a  few  to  consider  the  need  for  attaining  high  power  density.  Selection  of  penetrator  density,  velocity 
characteristics,  and  standoff  of  the  velocity  gradient  rods  are  now  available  for  parametric  evaluation. 

2.2  Linear  Velocity  Gradient  Rod 

The  case  of  the  linear  velocity  gradient  penetrators  presents  a  more  interesting  investigation.  The 
gradient  imposes  standoff  sensitivity  of  the  penetration  relationships.  Additionally,  the  high  velocity 
typically  resulting  from  shaped  charge  jets  is  seen  to  significantly  contribute  to  the  momentum,  energy 
and  power.  Less  obvious,  but  of  relevance  to  the  discussions  made  herein,  is  the  contribution  of  the 
material  densities.  The  density  ratio,  y,  is  found  to  contribute  to  the  spatial  terms  of  momentum,  Eq. 
(13),  to  the  2nd  power  and  energy,  Eq.  (14),  to  the  3rd  power.  An  illustration  for  the  energy  deposition 
is  given  in  Figure  4  for  a  linear  velocity  rod  having  mass  of  0.1kg  with  leading  tip  velocity  of  1 0km/s 
and  tail  velocity  of  3km/s.  Three  different  material  densities  are  selected  -  aluminum,  copper,  and 
tungsten  -  to  show  how  energy  is  deposited  according  to  y  dependency. 

The  low  density  materials  deposit  their  energy  in  the  initial  phase  while  the  high  density  materials 
carry  a  greater  percentage  of  their  kinetic  energy  deeper  into  the  target.  For  this  example,  the  A1 
deposits  over  2.5  times  the  energy  than  the  U  rod.  Recognize  that  this  is  only  an  input  quantity  and 
does  not  account  for  post  penetrator/target  interface  events  or  other  complex  target  response 
mechanisms.  However,  the  response  curves  are  useful  as  a  fundamental  criterion  or  screening 
parameter  for  selecting  a  penetrator  material  to  address  a  specific  target  material  and/or  penetration 
objective.  One  example  is  that  a  low  density  shaped  charge  liner  material  might  be  used  to  maximize 
the  entrance  hole  area  with  shallow  depth  for  mining  purposes. 
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Fig.  3.  Density  term  of  the  power  input,  Eq.  (5),  for  various  rod  materials  penetrating  a  steel  target. 
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Fig.  4.  Effect  of  penetrator  density  on  energy  deposition  for  a  linear  velocity  gradient  rod. 


(l+3y) 


D.E.  Lambert  /  International  Journal  of  Impact  Engineering 


2. 3  Energy  Deposition  of  a  Shaped  Charge  Jet 

The  effect  of  power  deposition  is  made  through  an  experiment  of  a  shaped  charge  jet  penetrating  a 
finite  thick,  but  severely  over-matched  piece  of  armor.  A  5-inch  diameter,  copper-lined  warhead 
(typically  used  for  very  deep  penetration)  was  fired  against  only  a  6-inch  thick  armor  plate  at  3 -charge 
diameter  (CD)  standoff.  A  radiograph  of  the  jet  at  150us  is  shown  in  Figure  5  with  the  tip  traveling  at 
9.54  km/s.  The  high-degree  of  overmatch  in  target  thickness  versus  its  total  penetration  capability 
resulted  in  a  violent,  explosion-type,  of  target  spallation  as  seen  in  Figure  6.  Following  the  1-D 
analysis  presented,  Eq.  (9),  and  the  virtual  origin  approximations  for  this  shaped  charge  design,  the 
velocity  of  the  eroded  tip  at  the  time  it  reaches  the  rear  surface  of  the  6-inch  target  is  8.51  km/s.  This 
segment  of  jet  length  has  then  proceeded  through  the  6-inches  of  steel  with  an  input  power  deposition 
of  8.5  MW. 

Photographs  of  the  target  plate  from  this  hypervelocity  penetration  experiment  are  given  in  Figure 
6.  The  rear  surface  has  a  noticeable  square-ended  cavity  with  the  surface  dominated  by  ductile  dimple 
and  pure  tension  features.  The  rapid  deposition  of  energy  and  the  intense  reflected  shock  from 
hypervelocity  impact  produced  such  severe  and  localized  spallation  event. 


Direction  of  travel  ^ 

Vtip=9.54km/s 


Fig.  5.  Radiograph  of  a  high  velocity,  copper  liner  shaped  charge  jet. 


Fig.  6.  Target  of  significant  over-matched  penetration  event.  Cross-section  cut  (far  right)  shows  severity  of  spallation  with 

power  input  of  approx.  8.5MW. 
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2.4  Hole-Volume  and  Energy’  Relationship 

Additional  analytic  analysis  is  possible  for  estimating  hole  volume  and  penetration  cavity  profiles 
using  the  work  of  Murphy,  et  al.  [6,7].  The  references  provide  a  basis  for  quantifying  the  energy  per 
unit  volume  of  target  removed  for  a  variety  of  materials.  The  volume,  V,  of  target  material  removed 
was  found  to  be  directly  proportional  to  the  jet  energy,  V  =  Elc,  with  c  being  a  constant  in  the  most 
basic  application  of  the  model.  The  observation  that  penetrator  material  was  only  a  second  order  effect 
on  hole  volume  reinforces  the  utility  of  this  paper’s  1-D  analysis.  More  complicated  forms,  including 
that  of  Perez  [8]  and  Szendrei  [9]  can  be  applied;  but  at  the  hypervelocity  regime  of  the  shaped  charge 
jet,  then  the  approximation  of  a  fixed  constant  is  suitable.  The  constant  proportionality  between  jet 
energy  deposited  and  target  volume  removed  is  seen  in  the  Figure  7,  extracted  from  Ref.  [6]. 


0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0 

Velocity  (kps) 


Fig.  7.  Energy-hole  volume  removed  from  Ref.  [6]  for  various  materials 


Murphy’s  relations  are  easily  incoiporated  into  the  differential  form  of  the  energy  relationships  to 
give  the  form  of  Eq.  (16).  Values  of  the  proportionality  constant,  c,  are  found  in  Figure  7  for  common 
metals  and  in  Table  1  for  some  geologic  materials  and  can  be  directly  applied  in  such  analysis. 


dV 

dp 


dE(d\ ^ 
dp 


dE 


with 


7  Ref  6 


d\ 

dE 


1 

c 


(16) 


Table  1 .  Energy  per  unit  volume  constant  for  various  geologic  materials  Rcf.[  10] 


Target 

c 

Material 

(J/cc) 

Concrete 

882 

Limestone 

1976 

Granite 

2409 
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3.  MATERIAL  BACKFLOW 


The  final  discussion  is  of  the  impact,  material  backflow,  and  interface  velocities  given  back  in 
Figure  1.  This  is  a  departure  from  the  energy-volume  relations  just  discussed,  but  is  relevant  to  the  1-D 
analysis  originally  posed.  Staying  in  the  coordinate  system  fixed  to  the  interface,  the  ratio  of  backflow 
material  (i.e.  material  turned  from  the  interface  impact)  to  the  impact  velocity  is  studied.  That  ratio  is, 
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backflow 


V. 


impact 


V  -2 u  _  y  - 1 

V  -u  / 


(17) 


This  ratio  delineates  the  penetrator/target  combinations  where  the  backflow  transitions  from  a  re¬ 
bounding  (opposite  direction  to  the  incoming  penetrator)  interaction  to  an  incoming  (same  direction  as 
the  penetrator).  Recall  that  this  is  with  respect  to  a  reference  system  attached  to  the  moving  interface. 
The  works  of  Allen  and  Rogers  [11]  describes  a  ‘secondary  penetration’  that  may  occur  and  the 
relationship  of  Eq.  (17)  and  the  following  discussion  compliment  their  hypothesis. 

The  response  curves  of  the  backflow-impact  velocity  ratio  are  shown  in  Fig.  8  for  penetrator 
material  densities  impacting  various  target  densities.  The  transition  points  (i.e.  material  is  leaving  the 
interface  at  the  same  rate  as  the  incoming  penetrator)  are  at  the  expected  values  where  the  impactor  and 
target  are  of  identical  density.  The  response  curves  also  support  the  phenomenon  commonly  observed 
in  experiments  having  high  density,  ductile  metals  (e.g.  tantalum)  penetrating  lower  density  targets  (e.g. 
A1  or  steel)  where,  for  specific  velocities  and  sufficient  target  depth,  the  penetrator  is  recovered  as  a 
‘thimble’  or  inverted  cylinder,  yet  totally  intact  as  in  Figure  9.  The  negative-valued  domain  of  Figure  8 
has  the  backflow  material  moving  in  the  same  direction,  albeit  lower  in  magnitude,  as  the  incoming 
penetrator.  The  unique  condition  of  post-recovery  of  a  fully  inverted  cylinder  also  requires  material 
ductility  and  flow  stress  to  turn  along  the  streamline,  as  well  as  have  sufficient  strength  to  resist  the 
strain  and  retard  the  gradient  along  the  inverted  penetrator. 


Fig.  8.  Backflow  velocity  response  curves 
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Fig.  9.  Diagram  of  the  inverting  cylinder  during  penetration 


The  sign  of  the  ratio  also  dictates  the  momentum  and  energy  exchange  to  the  penetration  interface. 
This  was  quantified  in  the  earlier  presentation  of  spatial  energy  and  momentum  deposition  terms  and 
graphs.  Penetrators  of  low  density  (with  respect  to  the  target  density)  have  greater  efficiency  on  the 
momentum  and  energy  exchange  with  an  analogy  made  to  perfectly  elastic  and  perfectly  plastic  impacts 
in  rigid  body  dynamics  studies. 


4.  CONCLUSIONS 

In  this  paper,  a  one-dimensional  analysis  of  fundamental  penetration  was  made  with  relationships 
that  characterize  the  energy  and  momentum  input  as  a  function  of  material  densities.  These 
relationships  are  useful  in  the  understanding  of  the  effect  of  material  density  and  the  ability  to  tailor  or 
dictate  desired  penetration  goals.  Low  density  penetrators  deposit  a  greater  percentage  of  their  energy 
in  the  initial  phase,  while  the  higher  density  materials  carry  their  energy  deeper  within  the  target.  An 
example  a  shaped  charge  jet  penetration  of  a  severely  over-matched  target  was  used  to  illustrate  how 
control  of  energy  deposition  might  be  exploited.  By  matching  penetrator  density  to  target  density,  the 
momentum,  energy  and  power  relationships  can  be  used  to  guide  first-order  objectives. 

The  analysis  is  ideally  suited  for  applying  energy  versus  hole  volume  data  and  models  from  the 
research  database.  Estimates  for  hole  profiles  can  be  directly  obtained  knowing  initial  characteristics  of 
rod  and  shaped  charge  jet  penetrators  along  with  the  analytic  equations  presented.  Material  strength 
and  complex  penetrator/target  interactions  were  not  addressed,  herein,  but  should  be  considered  for 
regimes  outside  of  the  hydrodynamic,  incompressible,  and  inviscid  domain  of  this  derivation. 

Further  discussions  were  made  on  the  interfacial  velocity  and  the  relative  nature  of  the  reverse 
streamline,  or  backflow.  The  relationship  should  be  considered  in  the  fundamental  understanding  of  the 
penetration  process  and  how  one  might  be  able  to  exploit  this  process  through  material  selection.  These 
relations  are  a  useful  addition  to  the  abundance  of  subject  matter  that  exists  relating  to  high  and 
hypervelocity  penetration. 
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